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The direct comparison of two different continuum discretization methods towards the solution of 
a composite particle scattering off a nucleus is presented. The first approach - the Continumm- 
Discretized Coupled Channel method - is based on the differential equation formalism, while the 
second one - the Wave-Packet Continuum Discretization method - uses the integral equation formu- 
lation for the composite-particle scattering problem. As benchmark calculations we have chosen the 
deuteron off a 58 Ni target scattering (as a realistic illustrative example) at three different incident 
energies: high, middle and low. Clear non-vanishing effects of closed inelastic channels at small and 
intermediate energies are established. The elastic cross sections found in both approaches are very 
close to each other for all three considered energies. 
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I. INTRODUCTION 

Historically the main progress in studying the elas- 
tic and inelastic (breakup) composite projectile scatter- 
ing off a heavy target was done within the Continuum- 
Discretized Coupled-Channel (CDCC) method [1-17] 
in the framework of the Schroedinger coupled-channel 
scheme. In the CDCC approach the total three-body 
scattering wave function (in initial versions of the method 
the projectile was assumed to be a two-fragment nucleus) 
is expanded in a complete basis of the fragments rel- 
ative motion of the projectile and the unknown three- 
dimensional elastic and breakup channel wave-functions 
of the projectile cm. motion are obtained by solving a set 
coupled-channel equations. This approach has been suc- 
cessfully applied in the field of nuclear reactions to find 
the scattering amplitudes of many particular processes 
for various types of projectiles, like deuterons, 6 Li, 6 He, 
7 Li, 12 C etc. However, a strong discussion about the va- 
lidity of the method arose [H, [H], HI] because there were 
some unclear issues related to the treatment of the to- 
tal wave function breakup components in the asymptotic 
region. So, to check the accuracy and reliability of the 
results attained within this approach, the convergence 
of the elastic and breakup amplitudes (with variation of 
the momentum bin widths and their number) was stud- 
ied in detail for various approximation schemes [H . This 
study has demonstrated clearly that the approach leads 
to convergent results for any reasonable choice of the dis- 
cretization scheme. In particular, it was shown in Ref. [9( 
(for the case of deuteron scattering off a 58 Ni target as a 
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test case) that the employment of a finite number of exact 
two-body scattering wave functions as a basis for the in- 
ternal projectile motion with energies fixed at midpoints 
of the bins (the midpoint method) leads unexpectedly in 
general to slower convergence for the elastic and breakup 
amplitudes than the usage of the two-body wave func- 
tions averaged inside energy bins (the average method). 
This result seems rather unnatural because the use of av- 
eraged wave functions (i.e. wave-packets) as a basis for 
the internal projectile motion leads to a wrong asymp- 
totic behavior of the coupling potentials and scattering 
wave functions in breakup channels, while the employ- 
ment of the exact wave functions of the projectile internal 
subHamiltonian seems from the first glance to be more 
adequate. 

Thus, the energy averaging method for the projectile 
internal continuum wave-functions seems to assist to get 
faster convergence in the CDCC calculations. The aver- 
aging leads to normalized (i.e. of L2-type) wave-functions 
(j>i(r), depending on the relative coordinate between the 
two fragments. At the same time, the motion along the 
second, projectile center of mass, coordinate R is consid- 
ered in the CDCC approach without any averaging on 
energy. So, the three-body wave function for the i-th 
excited channel has an outgoing asymptotic behavior 

*(r,R) ► ^(r)Sf^. 

r<HR 

Therefore above the three-body breakup threshold the 
total wave function ^(r, R) vanishes at r — ► oo (simi- 
larly to the bound-state function) while the dependence 
of the total scattering wave function upon the projectile 
center of mass coordinate has the pure outgoing wave 
asymptotic behavior. Such an asymmetry for the motions 
along the two coordinates at ener gies above the breakup 
threshold looks not fully justified [l6l \vA HH because in 
this energy region the proper asymptotic behavior over 
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all spatial variables might be important. 

However, in case of the scattering of a weakly bound 
projectile (e.g. deuteron or 7 Li) at relatively high incident 
energies (E > 40 MeV) the motion along the projectile 
cm. coordinate R proceeds much more faster than the 
motion along the relative coordinate r (if the projectile 
cm. energy is much larger than its binding energy). In 
this case the breakup of the projectile takes place mostly 
outside the nuclear held of the target. In other words, 
the projectile is first excited into its internal-continuum 
states and then (already outside the reaction region) the 
excited projectile breaks up into its fragments. Within 
this picture, the breakup of the projectile in the tar- 
get field can be treated as a particular case of inelastic 
scattering into the projectile continuum states. Then a 
proper averaging over energy bins for the fragments mo- 
tion in the discretized continuum gives just a discrete 
representation for an exact spectral density function of 
this subsystem. As a result, this discretized density func- 
tion enters all the elastic and breakup amplitudes. Thus, 
the applicability condition for the treatment of the above 
processes within the CDCC method should be as fol- 
lows: £totai 3> ^binding- However when the projectile 
kinetic energy Btotal is comparable or even a bit larger 
than its inner binding energy the other asymptotic chan- 
nels should be taken into consideration and thus, the 
coupled-channel Schroedinger equation formalism should 
be replaced by the respective Faddeev approach. 

Thus, within its scope of validity, the CDCC approach 
can be considered as a rather useful and practical working 
method to handle with the realistic problems of compos- 
ite projectile scattering off a target. However it would be 
very important to test carefully the accuracy and predic- 
tive power of the CDCC results in some fully independent 
way based on a different approach. The present paper is 
devoted just to a careful comparison between the CDCC 
and one alternative method for some typical problem in 
the field. 

The alternative approach which we have chosen for our 
study is a new integral equation approach based on a to- 
tal continuum discretization [Til, 120 . l2ll . |22T ] which has 
been developed by two of the present authors some time 
ago. The approach is based on the Wave-Packet Con- 
tinuum Discretization (WPCD) technique which allows 
to convert the three-body integral equations into the re- 
spective matrix equation with the matrix elements easily 
calculable in an analytical way. In the WPCD method 
the motion along both radial coordinates r and R is con- 
sidered symmetrically, i.e. the respective averaged wave 
functions (stationary wave packets) vanishing at infinity 
along both coordinates are used. Such a fully discretized 
integral approach can be used because in the expression 
for the elastic scattering amplitude the total three-body 
resolvent operator G(E) is bracketed from the left and 
from the right by fast decaying factors over all the spa- 
tial coordinates: 

A cl (E) = {^^ { ~\E-e )\V+VG(E)V\^^£\E-e )), 



where \<f>o,<ipQ (E — e )) is the initial (final) projectile 
three-body wave function, so that the "external" inter- 
action, V, decaying along the cm. coordinate R and 
the projectile bound-state wave functions 0o( r ) decay- 
ing along the internal coordinate r cuts effectively the 
asymptotic parts of the total resolvent G{E). Thus, in- 
stead of the exact three-body resolvent operator G(E), 
its respective wave-packet finite-dimensional approxima- 
tion G(E) can be favorably used [2(J 22]. Therefore, in 
this integral wave-packet approach we do not need any 
asymptotic parts of scattering wave functions. All that is 
required is just the inner parts of the wave functions but 
with the proper normalization which can be well approx- 
imated via the stationary wave packets in the WPCD 
approach. It should be stressed that the effective cutoff 
of peripheric three-body resolvent parts does not depend 
on the total energy. Hence, if to neglect the stripping 
channels (similarly to the conventional CDCC method) 
and to treat the elastic and breakup channels only, the 
wave-packet method looks fully justified at low and high 
energies as well. 

Such a 'totally discretized' integral technique has sev- 
eral evident advantages as compared to other formula- 
tions of the same problem. In addition to a convenient 
matrix scheme for a calculation of both the elastic and 
breakup amplitudes, it also allows to construct in an ex- 
plicit form an optical (energy-dependent and non-local) 
effective potential for the projectile-target interaction 
which accounts for the intermediate inelastic processes 
properly [2p| . [22| . Another distinctive feature of the inte- 
gral approach is a proper and straightforward inclusion 
of closed channels which are treated in this approach on 
the same footing as the open channels (see especially the 
section III). 

So, to make the comparison between two approaches 
maximally broad we considered the deuteron elastic scat- 
tering by 58 Ni at three different incident energies: high, 
intermediate and low, where closed channel contribu- 
tions are expected to be quite different. This reaction 
has been also used to compare the CDCC and Faddeev 
approaches in a recent work [23j . showing very good 
agreement between both methods. It will be vividly 
demonstrated that the contribution of the closed chan- 
nels, which are commonly neglected in CDCC calcula- 
tions P, [3, S H S 0, 0, H0, raises essentially 
when the projectile energy gets low. In principle, closed 
channels can be incorporated within the CDCC scheme, 
however, in this case the solution of the coupled equa- 
tions becomes numerically unstable. In this situation, 
it is more convenient to solve the problem using the R- 
matrix method [24j . This is an equivalent way of solving 
the coupled equations, but with the advantage of being 
more stable numerically as compared to the traditional 
radial stepping methods. However, this reformulation 
demands more computational resources than the con- 
ventional coupled-channel scheme [ill ]. Thus, we made 
a comprehensive comparison for both open and closed 
channel contributions between the R-matrix CDCC and 
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WPCD approaches at different energies. 

The paper is organized as follows. In Section II a brief 
description of the CDCC and WPCD approaches is pre- 
sented. Section III contains numerical results obtained 
with these two approaches and the discussion is placed 
in Section IV. Finally, the summary of the paper is given 
in Section V. 



II. BRIEF DESCRIPTION OF THE CDCC AND 
WPCD APPROACHES 

The Hamiltonian describing the scattering of the two- 
fragment composite particle {be} off the target nucleus 
A is taken in the form: 

H = h bc (r)+hc(R)+V b - A (E la . h /2)+V c - A (E lab /2)+AV c , 

(1) 

where hb c is the subHamiltonian of the internal b — c 
motion (acting on the fragment relative coordinate r), 
he — the subHamiltonian of the projectile center of 
mass asymptotic motion including long-range point-like 
Coulomb interaction (acting along the cm. coordinate 
R) , V c -a and V b -A are energy dependent optical poten- 
tials for the fragment-target system taken at half the 
incident particle energy E\ a \j (for the nearly equal mass 
fragments). Finally, AVc is the additional short-range 
Coulomb interaction (acting along the cm. coordinate) 
which is caused by the finite charge radius of the target 
nucleus (due to charge distribution in the target nucleus). 
To obtain the elastic scattering and breakup amplitudes 
for the above scattering problem with the three-body 
Hamiltonian ([T]) one has to solve either the three-body 
Schroedinger equation with proper asymptotic bound- 
ary conditions or the system of the respective three-body 
Lippmann-Schwinger equations. 



A. Continuum discretized coupled channel method 
(a differential equation formulation). 

In the CDCC approach the continuous spectrum of 
the hbc subHamiltonian is discretized by dividing the 
continuous momentum distribution into a finite num- 
ber of non-overlapping bins ki]fL 1 with correspond- 
ing averaged (within each bin) continuum wave functions 
{\<t>i)}fLi and respective averaged 'channel eigenenergy' 
values e* = (4>i\hb c \4>i} . To simplify the notation we omit 
here the angular variables but they are assumed to be in- 
cluded. Then, the total three-body wave function of the 
Hamiltonian H is expanded over the set of averaged wave 
functions describing the projectile internal spectrum (in- 
cluding the projectile bound state |(/>o) with the binding 
energy 6q) : 

N 

\^(E))=J2\^Xi), (2) 

i=0 



where (R|xi) = Xi(R) are the channel wave functions 
which define the elastic (i — 0) and breakup (i ^ 0) am- 
plitudes. The bin wave functions {\(j)i}}iLi are typically 
constructed by averaging the true continuum states, </>(fc) 
within the bin interval: 

|&) = -L t /(fcMfc))dfc, i=l,...,N (3) 
V A » Jki-i 

where Aj is a normalization constant and f(k) is a 
weight function. For non-resonant continuum the com- 
mon choice is f(k) = 1, in which case Aj = ki — fcj_i. 
Applying the expansion @ to the initial Schroedinger 
equation for the Hamiltonian |T]) one gets a system of 
coupled differential equations for the unknown functions 
Xi(R): 

N 

[/ic(R)+AVb(R)+V«(R)-(£-ei)]x<(R) = f«'Xi'(R), 

i'=0 

(4) 

where V«'(R) = (<pi\Vb-A + V c -A\<k') 1 are coupling po- 
tentials. This system is solved with the following bound- 
ary conditions [9(: 

X i(R)~I(Ki,R)5 i0 - x l*ks ifi O{Ki,R), (5) 

V AO 

where I and O are the Coulomb incoming and outgoing 
waves, Ki is the cm. wave number which is related to 
the 'channel eigenenergy' e* by energy conservation, i.e. 

E = e* + ^2m" ( nere M is deuteron-target reduced mass) 
and Sifi are S-matrix elements for the elastic (i = 0) and 
inelastic (i ^ 0) scattering 2 . 



B. Wave-packet continuum discretization method 

In the WPCD approach one applies an analogous dis- 
cretization of the continuous spectrum of the h bc sub- 
Hamiltonian but making use of energy bins: [e,-_i, £i\f =1 3 ■ 
In a similar way, one defines a set of stationary wave 
packets (WP) as integrals of the exact scattering wave 



1 The integral is taken over the relative r coordinate only. 

2 For closed channels (ej > E) the boundary condition should be 
defined with an exponentially decaying and rising functions. In 
the standard CDCC approach the inclusion of closed channels 
causes numerical instabilities due to the unavoidable appearance 
of exponentially rising asymptotic parts, so that the closed chan- 
nels are usually not taken into account in these calculations. To 
avoid this problem the combination of the CDCC and R-matrix 
approaches can be used. 

3 It is possible to make just the same energy and momentum bins 
in the WPCD and CDCC approaches using the interrelation 

6j = 2 * , where m is the reduced mass of the two projectile 
constituents. 
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functions \4>{E)) of the h bc subHamiltonian over the re- 
spective energy bins: 




\<t>(E))dE, i = l,...,N (6) 



where D t = e, — ej_i are the bin widths and the cor- 
responding 'channel eigenenergies' arc just the bin mid- 
points: e* = |(e-i-i + £j)- The set of these WP states 
should be supplemented with the bound state \Zq) = 
\4>o). The main difference in this discretization proce- 
dure with respect to the conventional CDCC one is that 
we use here a full set of WP states which includes both 
open and closed channels. 

Besides that we build a partition of the continuous en- 
ergy spectrum of the Coulomb subHamiltonian he act- 
ing along the projectile cm. coordinate R and employ 
the corresponding set of Coulomb stationary wave pack- 
ets (CPs) constructed by averaging the regular Coulomb 
continuum wave functions \ip c (E)) over energy bins 
[£j—i,£j]j = i_: 

\Xf) = -L= t \^ c (E))dE, (7) 

where Aj = £j — £j-i are the bin widths. Finally, we 
construct the three-body wave-packet (TWP) basis set 
for the elastic channel asymptotic Hamiltonian, 

H bc = h bc ®h c , (8) 

as the direct product of the states \Zj) and \Xj) (includ- 
ing all the required partial wave couplings): 

\Sij) = \Zi,Xj), i = 0,...,N; j = l,..., M, (9) 

Such a three-body packet basis (which is also of the L2 
type) is very convenient for the subsequent investigations. 
One of its main advantages is that the matrix of the chan- 
nel resolvent G bc = [E + iO — H^]" 1 is diagonal in this 
basis. Moreover, the corresponding matrix elements of 
the three-body resolvent operator have an explicit ana- 
lytical form [20I |22| and they depend on the spectrum 
discretization parameters (i.e. and £ j values) only. 

Due to the fact that the interactions between the pro- 
jectile constituents and the target are given by complex 
optical-model potentials the contribution of the inter- 
mediate rearrangement (viz. stripping) channels to the 
elastic scattering or breakup should be of minor impor- 
tance 4 . Therefore, in order to treat the projectile elastic 
scattering in the leading order, one can take into account 
only those intermediate breakup states whose asymptotic 
behavior are defined by the three-body channel Hamil- 
tonian H bc 5 . So, the total three-body scattering wave 



function \^(E)) in this case can be found from a single 
Lippmann-Schwinger (LS) equation: 

|tf (£)) = \* d (E)) + G bc (E)V bc \V(E)), (10) 

where V bc = V b -A + V c -a + &Vc is the "external" short- 
range interaction and \^ d (E)) = \Z ,ip c (E — eg)) is the 
outgoing wave function for the projectile in its ground 
state. The long-range Coulomb part for the point-like 
Coulomb interaction between the projectile and target 
is included into the G bc operator. Next, we make the 
projection of Eq. (fT0|) onto the wave-packet subspace and 
expand the total wave function in the three-body wave- 
packet basis for the Hamiltonian H bc : 

mE))n\*(E))=J2Ci j \S ij ). (11) 

ij 

Then Eq. (|10|) can be rewritten in a finite-dimensional 
form 

\$(E)) = \S 0jo ) + G bc (E)V\$(E)), (12) 

where \Soj ) = \Zo,X!?) is the packet state correspond- 
ing to the outgoing wave function \&d} (the index jo is 
defined by the energy rule: E e [£j -i + e 0'^jo + e o])- 
The interaction operator V with a hat is the wave packet 
projection of the external interaction V bc 

V^Y. {Si'j'KSijlVw, (13) 

whose matrix elements in the TWP basis can be inter- 
preted as generalized coupling potentials: 

Vwj, = (S^V^a + Vc-a + AVb|Si'iO (14) 

and may be calculated either analytically or numerically. 
Thus, in this approach solving the three-body LS equa- 
tion is reduced to the solution of the simple matrix equa- 
tion (fT2"l) 6 . 

The composite-particle elastic scattering amplitude 
can now be written formally as the following diagonal 
(on-shell) matrix element: 

A cl « [( V ~ 1 -y- 1 ]o J -o,o J -o ) (is) 

where matrices of the respective operators in three-body 
WP basis are denoted by the bold letters, while A J0 is 
the bin width of the he continuum partition. Moreover, 



4 Due to the imaginary potentials there are no bound states in 
the rearrangement channels, and also the kinetic energy of the 
projectile is assumed to be relatively high. 

5 This is just the same assumption as in the CDCC approach. 



6 It should be stressed that the matrix equation 11121 in our case 
is strongly distinguished from that which can be deduced from 
a direct matrix reduction of the three-body LS equation using 
some quadrature mesh-points because all the matrix elements in 
our case have been averaged over the energy bins corresponding 
to both the active coordinates. 
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the breakup amplitude can be readily found as a non- 
diagonal (off-shell) matrix element of the same matrix 
[V" 1 - Gbc]- 1 in the WP basis [13, HI- The above ma- 
trix is essentially a finite-dimensional analog of the exact 
transition operator matrix. 

Here it is appropriate to summarize the main distinc- 
tive features which distinguish the wave-packet approach 
from the conventional CDCC method: 

(i) The wave-packet basis for the projectile bound and 
continuum internal states is chosen to go up to high 
excitation energies, so that not only open-channel 
but also closed-channel contributions are fully in- 
cluded into the WP three-body calculations. 

(ii) The 'internal' and 'external' motions of the pro- 
jectile are taken into consideration in the WP tech- 
nique on 'equal footing', i.e. using the bases of wave 
packets averaged over the respective energy bins for 
the motions along both independent radial coordi- 
nates. So that the total three-body resolvent in the 
transition matrix elements is cut effectively at large 
distances by the interaction potentials and also by 
the bound-state functions of the projectile in the 
initial and final states. 

(iii) The above mentioned cutoff of the total three-body 
resolvent is intimately related to the integral for- 
mulation of the WP approach, in which a proper 
asymptotic behavior of coupling potentials and also 
an explicit matching procedure between inner and 
asymptotic wave functions are not required. 

(iv) The kernel of the respective matrix equation in the 
WP approach includes rather smoothed (on energy) 
matrix elements due to the integration over energy 
bins for the motions along all active coordinates. 
This important feature of WP matrix kernels can 
readily be extended to the solution of the general 
Faddeev equations using the WP-technique [2(| . 

(v) The integral formulation of the three- and few-body 
scattering problem makes it possible to reformulate 
easily the initial problem in terms of the Feshbach 
projection operator technique (22J. This reformu- 
lation allows to reduce the elastic scattering of a 
composite projectile to a simple potential scatter- 
ing by a non-local and energy-dependent potential 
and moreover to construct this potential in an ex- 
plicit analytical form. 



III. COMPARISON BETWEEN THE CDCC 

AND WPCD: THE BENCHMARK 
CALCULATIONS AT DIFFERENT INCIDENT 
ENERGIES 

As a general numerical test of the two methods we 
have chosen the popular case of the elastic scattering of 
the deuteron off 58 Ni target at three different incident 
energies, namely, E\ a h — 80, 21.6 and 12 MeV. These 
detailed calculations allow not only to make more wide 
comparison between two approaches but also to deduce 
from this comparison the role of the closed channels at 
different energies. 

For the sake of simplicity, the n-p interaction is as- 
sumed here to have the Gaussian form 0]: 

V(r) = -V e- pr \ Vo = 66.99 MeV, /3 = 0.415 fnT 2 . 

(16) 

Also, keeping in mind the comparison purposes between 
the two alternative methods, we restrict ourselves here 
to s- waves for the n-p relative motion 7 . 

The nucleon-nucleus interactions are represented by 
local optical-model potentials of Woods-Saxon type. 
In particular, we employ here the potentials from the 
Koning-Delaroche global fits [25| evaluated at half of the 
deuteron incident energy. The mass values m n — 1.0087 
amu (for the neutron), m p = 1.0078 amu (for the proton) 
and mNi = 57.9353 amu (for the target) are adopted here. 

A. Choice for parameter values of the momentum 

distribution in the CDCC method 

Following the standard procedure, in the CDCC calcu- 
lations the continuum was truncated by defining a max- 
imum excitation E max , and divided into ./V energy bins, 
evenly spaced in the linear momentum. For each bin, 
a representative wave function is constructed according 
to Eq. with f(k) = 1. The number of bins N was 
determined in order to achieve convergence of the elas- 
tic angular distribution. For the two lower energies, 
Slab = 12 MeV and E lah = 21.6 MeV, the calculations 
include both open and closed channels. Thus, the con- 
tinuum was represented by a set of N opon bins up the 
maximum available cm. energy, and by N c \ ose d bins from 
this energy to -E max , and hence N = N opcn + N c i oscd . 

B. Choice for parameter values of the wave- packet 

bases for the WPCD calculations 

To simplify the construction of the n-p relative mo- 
tion WPs \Zi) and to get simultaneously the bound state 



So, it would be very instructive to make a direct com- 
parison between the WPCD and the conventional CDCC 
results for some standard test problem. 



7 Higher n-p partial waves would be certainly important to com- 
pare the theoretical results with the respective experimental 
data. 
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FIG. 1: Energy (a) and momentum fcj (b) distributions for 
the n-p continuum discretization constructed with a Gaussian 
basis on the generalized Tchebyshev grid (|18p . 



of the exact Coulomb packets of the dimension M con- 
structed also on the Tchebyshev grid with t = 1, which 
corresponds to the energy distribution: 

Ej = E Q tan r^' 1 ^ , j = l,...,M, (19) 

where £j are the end points of the energy bins (£q = 0) 
and Eq is the incident energy of the projectile in the cm. 
frame. 



C. Results of calculations 

In this section we will present the particular results for 
the differential elastic cross sections calculated at differ- 
ent incident energies. 

1. £ lab =80 MeV 



wave function \Z ) we apply here a single diagonaliza- 
tion procedure for the hb c subHamiltonian on the simple 
Gaussian basis \i/> n )n=i> i- e - we assume that the n-p wave 
packets are expanded as follows: 



\Zi) 



N G 

E 



A n \i> n ), ip n (r) = Bj, 



(17) 



where B n are normalization factors. The scale parame- 
ters a n are defined on the generalized Tchebyshev grid 



a 



tan 



n{2n- 1) 
4N n 



ii = 1, . 



(18) 



where ao is the common scale parameter and t defines 
the distribution sparseness for the Gaussian 'frequencies' 
a n . As t become larger the mesh points a n are simul- 
taneously extended to lower and higher values, so that 
more longer-range and shorter-range Gaussian compo- 
nents appear in the basis. The energy distribution e* 



(and the corresponding momenta ki 



2me* 



obtained 



with such a basis occurred to be highly non-homogeneous 
(i.e. non-equidistant), but it covers practically all the en- 
ergy spectrum and provides well converged results for 
the scattering amplitudes [l9[ . Figure [1] shows the en- 
ergy and corresponding momentum distributions for the 
following parameter values chosen for the grid: Ng = 20, 
t = 2, ao = 0.1 fm~ 2 . In practical calculations we found 
that it is quite sufficient to include only the first N basis 
functions with maximum bin energy ejv less than some 
chosen maximal energy value E max of the included closed 
channels. For example, it is clear from Fig. Q]that a ba- 
sis containing the first N = 13 functions is enough for 
£max = 100 MeV. 

As for the WP basis to describe the n-p pair cm. mo- 
tion we employ (for all required partial waves L) a set 



At this relatively high energy, we found that the 
breakup effect contribution is rather small, producing 
just a minor correction to the folding model result (i.e. 
the calculation omitting continuum channels). Also it 
was shown previously that only about half of the open 
channels (i.e. E max w 40 MeV) play a significant role at 
this relatively high energy. So one can conclude that the 
effect of closed channels can be ignored here. In Fig. [2] it 
is shown the comparison between the elastic differential 
cross sections obtained within the CDCC, the WPCD 
and the folding model. For these calculations the WP 
bases with dimensions N — 13 (for n-p relative motion) 8 
and M — 500 (for deuteron cm. motion) are used. The 
number of partial waves required for the deuteron-target 
relative motion is L max = 45. For the CDCC calcula- 
tions, the continuum was truncated at -E max = 70 MeV 
and divided into N — 20 bins, uniformly distributed in 
the linear momentum k. The set of coupled equations 
were solved using the Numerov method, and matched to 
their asymptotic solution at 20 fm. These calculations 
were performed with the code FRESCO [TT[ | . 

As it follows from Fig. [2] the results of both discretiza- 
tion methods are in very good agreement for all angles. 
It should be emphasized that despite the highly non- 
homogeneous bin distributions for the internal n — p sub- 
system as well as for the deuteron-target subHamiltonian 
lead to well converged results of the WPCD calculations. 
So, it is not necessary to use in practical calculations just 
equidistant momentum bin distributions as it is usually 
done 0, 0, [H • Moreover, the application of some special 
energy distribution allows to decrease significantly the 
number of the packet basis functions required for conver- 
gence. Thus, to summarize, despite of the quite different 



The Gaussian basis 03 with N G = 20, t = 2, ao = 0.1 fm -1 
has been used for this calculation. 
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FIG. 2: Comparison for the differential elastic cross sections 
of the deuteron off 58 Ni-target scattering (plotted as a ra- 
tio ratio to the Rutherford cross section ctr) at i?i a b = 80 
MeV calculated with the CDCC (dotted line) and WPCD 
approaches (full line). The dashed curve corresponds to the 
single folding model calculation. 



bin distributions considered in the CDCC and WPCD 
approaches and rather different formalisms both meth- 
ods provide results in nice agreement to each other at 
E lah = 80 MeV. 



2. -Biab =21.6 MeV. The effect of closed channels 

At this intermediate energy, the influence of closed 
channels (CC) on the elastic scattering amplitude be- 
comes visible. Our calculations have shown that the main 
part of this effect is seen at backward angles. Though the 
CC effect has been commonly neglected in the conven- 
tional CDCC calculations [1 1, MM !, !, 0, S 0, 03 it 
is worth to study carefully the convergence of the elastic 
cross section when the number of closed channels and the 
corresponding maximum energy £ max are increased. 

To prepare this calculations in the WPCD scheme we 
have used the three-body WP sets of dimensions N = 15 
(constructed with the Gaussian basis (jTTJ) with Nq = 30, 
t = 3, ao = 0.1 fan -1 ) and M = 150. The number of par- 
tial waves required at this energy is L max — 25. To in- 
corporate the closed channels into the CDCC framework 
a special combination of the coupled-channel scheme and 
the R-matrix approach has been used [24] . In contrast to 
the Numerov method, where the radial functions Xj(R-) 
are obtained by direct integration of the coupled differ- 
ential equations (U]), in the R-matrix method a basis set 
of energy eigenstates is first obtained by solving the diag- 
onal parts of the coupled equations, with the basis func- 
tions all having fixed logarithmic derivatives at a given 
distance R m . Next, the unknown functions Xi(R) are 
expanded in this basis. Although this procedure is com- 
putationally more demanding, it leads to more stable re- 



sults when closed channels are taken into account. In 
order to achieve convergence at this incident energy, we 
required N open ~ 20 and 7V c i ose< j ~ 15 bins, for a maxi- 
mum excitation energy .E max ~ 140 MeV. The R-matrix 
radius was fixed to R m = 30 fm. 

In Fig. [3] the comparison between the R-matrix CDCC 
and the WPCD elastic cross sections calculated at almost 
equal E max values is presented. For E max = 20 MeV, 
there are only open-channels and the conventional Nu- 
merov procedure has been employed to solve the CDCC 
equations (see Fig. (3^). For the other values of E max 
(Fig. [3)3 and c), the R-matrix method has been employed, 
while the WPCD technique has not been changed. It is 
evident that the results for the differential cross sections 
found with the WPCD and R-matrix CDCC methods 
with similar £ mM values are almost identical. So, it can 
be concluded that the details of the discretization distri- 
butions in both methods do not play an important role 
in the convergence and the most significant factor here is 
just the i^max value. 

The converged elastic cross sections (calculated with 
the highest value, E max = 139 MeV) are presented in 
Fig. [H The curves corresponding to both discretization 
methods are hardly distinguishable. In this figure we 
have included also the calculation corresponding to the 
folding model calculation (dashed line) . It is evident that 
breakup effects are large at this incident energy. 

In Fig. [5] we plot the elastic ^-matrix elements obtained 
with both discretization methods. It can be seen that 
they are in very good agreement for all partial waves. 

From the present results it becomes clear that the CC 
effect must be taken into account at these intermediate 
energies. Furthermore, the required maximum energy of 
the CCs is rather high, viz. E max ~ 100 MeV. 



3. E lBb =12MeV 

We finally present in this section the calculations at 
the relatively low energy of 12 MeV. In this case the 
dimensions of the WP sets N = 15 (constructed from the 
Gaussian basis (fT7|) with Nq = 30, t = 3, ao = 0.1 fm -2 ) 
and M ~ 150 are sufficient, and the maximum deuteron 
cm. angular momentum value is only L max = 15. For the 
CDCC calculations we used N open = 15 and N c i osc( i = 10 
(for E max ~ 80 MeV). Again, the R-matrix technique 
was used to solve the coupled equations, with a matching 
radius of R m = 30 fm. 

The elastic cross sections calculated with the WPCD 
and the R-matrix CDCC approaches for different E maK 
values are shown in Figs. [S^, b and c. The compari- 
son between the converged (with respect to E max values) 
cross sections calculated with both methods along with 
the folding model results are presented in Fig. [7J 

It is evident that the influence of the CC is the largest 
here and this effect is seen not only at backward angles 
but also at intermediate angles (0 c . m . ~ 90° — 100°). The 
small discrepancies between the WPCD and R-matrix 
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FIG. 3: Elastic d+ 58 Ni cross section at £i ab = 21.6 MeV cal- 
culated within the CDCC (dotted line) and WPCD (dashed 
line) approaches for different maximum values £ ma j< of the 
excitation energy: a) _E max = 20 MeV, b) _E max = 38 MeV, c) 
-Bmax = 80 MeV. The full lines correspond to the converged 
CDCC and WPCD results, which are almost indistinguishable 
in these plots. 

CDCC results observed at most backward angles might 
be ascribed to the different n-p continuum discretization 
parameters, but we can not rule out that they can be re- 
lated to numerical inaccuracies of either of the two meth- 
ods. In any case, this effect, which has not been seen in 
the two previous cases, begins to play some role at rather 
low incident energies. 

IV. DISCUSSION 

The close agreement between the WPCD and R-matrix 
CDCC approaches found in the present work can be in- 
terpreted as a good indicator that the basic assumptions 
(or approximations) adopted in the WPCD method, in 
particular the employment of orthonormalized stationary 
wave packets instead of exact non-normalizable three- 
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FIG. 4: Converged results (with respect to the maximum en- 
ergy -E m ax) for the differential elastic cross section (relative 
to Rutherford) for d+ 58 Ni at E uh = 21.6 MeV, calculated 
with the CDCC (dotted curve) and WPCD (solid curve) ap- 
proaches. The dashed line corresponds to the folding calcula- 
tion, in which continuum channels are not taken into account. 
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FIG. 5: Argand plot for the elastic S-matrix elements for the 
reaction d+ 58 Ni at _Ei a b = 21.6 MeV calculated within the 
CDCC (open circles connected by dotted lines) and WPCD 
(black squares connected by solid lines) methods. 

body scattering wave functions are well justified and cor- 
rect. This fact implies some further consequences in the 
whole many-body scattering theory. Actually, the us- 
age of the total continuum discretization on the basis 
of normalizable (of L 2 type) stationary wave packets is 
equivalent to the formulation of the many-body scatter- 
ing problem in a finite multi-dimensional box and leads 
a treatment of the original scattering problem analogous 
to the much more simple bound-state problem. In turn, 
the last conclusion allows to treat the many-body scat- 
tering problem in similar way as the many-body bound 
states, e.g. by expanding the L2 _ type stationary wave 
packets (constructed from the exact unknown scattering 
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FIG. 6: Elastic d+ Ni cross section at Slab = 12 MeV cal- 
culated within the CDCC (dotted line) and WPCD (dashed 
line) approaches for different maximum values of the excita- 
tion energy, £ max : a) fi malI = 10.2 MeV, b) £? max = 24.5 
MeV, c) -Emax = 36 MeV. The full curves correspond to the 
converged CDCC results. 
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FIG. 7: The same as in Fig. H but for E lsh = 12 MeV. 



wave functions) on some convenient Li basis widely used 
in bound state problems (for example, the many-body 
harmonic oscillator basis). The first results which two 
of the present authors [201 ] have obtained with the appli- 
cation of this technique for the solution of the Faddeev 
equations above the three-body breakup threshold seems 
to support this expectation. 

Another nice feature of the WPCD approach is the 
possibility to reformulate in a straightforward way the 
whole problem in terms of the Feshbach projection tech- 
nique [22j | and to include all the breakup channels into the 
respective Feshbach polarization non-local potential. Af- 
ter this, one can readily use the full applicability of the 
WPCD approach to the resulted non-local interactions 
(in contrast to the conventional coupled-channel method) 
and to reduce in a standard way the complicated many- 
body scattering problem to the simple matrix equations 

it 

On the other hand, the very close agreement between 
the conventional CDCC approach and the fully alter- 
native wave-packet integral method should remove any 
doubts reported in literature [lfl [H| concerning the re- 
liability and convergence properties of the Schroedinger 
coupled-channel scheme in the framework of the three- 
body scattering problem. 

V. CONCLUSION 

In this work, we have compared two essentially differ- 
ent approaches based on the general idea of continuum 
discretization for solving three-body scattering problems. 
The CDCC method was proposed to solve the coupled- 
channel problem in the framework of the Schroedinger 
equation formalism. This approach uses a discretization 
of the projectile internal continuum only, along with an 
explicit matching between inner and asymptotic scatter- 
ing wave functions to extract the multichannel S'-matrix. 
On the other hand, the WPCD approach is based on the 
total three-body continuum discretization which allows 
to use the integral equation formalism of the scattering 
theory. The special wave-packet technique allows to con- 
struct finite-dimensional analogs for all the scattering op- 
erators and to use only inner parts of wave functions in 
order to find the observables. 

The direct comparison of both methods leads to the 
following conclusions: 

(i) The two methods of continuum discretization dis- 
cussed in this work produce very close results for 
the composite particle elastic scattering at rather 
high and also at low energies. 

(ii) Different types of discretization of the projectile 
internal subHamiltonian continuum used in both 
methods do not lead to visible discrepancies at high 
and intermediate energies. Small differences in the 
elastic cross sections arise only at low incident ener- 
gies. An employment of a highly non-homogeneous 
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dicretization distribution in the WPCD approach 
still leads to well-converged results. 

(iii) The effect of closed channels is not seen at high in- 
cident energies. This effect begins to play a visible 
role at intermediate energies (at the most backward 
angles) and can not be ignored as it was assumed 
in some previous CDCC calculations. The conver- 
gence of the calculated cross section when increas- 
ing the number of included CCs depends only on 
the maximum excitation energy and is not sensitive 
to the details of the continuum dicretization. 

(iv) The effect of CC on the elastic scattering becomes 
more significant when the incident energy of the 
projectile decreases. 

(v) The very close agreement between the results ob- 
tained with these two different approaches at var- 
ious energies, from low to high, seems to indicate 



that the three-body wave-packet technique can be 
considered as a rather reliable and numerically ac- 
curate method for calculations of composite projec- 
tile scattering by nuclei. Simultaneously, one can 
confirm again that the conventional (or R-matrix) 
CDCC method leads to quite reliable and well- 
converged results for the elastic and inelastic scat- 
tering of composite projectiles. 
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